function [] = fig_sec11_7_modularCapacitiesAndRouting_dualDecomp()

close all;

objFunc = load ('data/modularCapacitiesAndRoutingDualDecomp_objFunc.txt');
bestObjFunc = load ('data/modularCapacitiesAndRoutingDualDecomp_bestObjFunc.txt');
lowerBound = load ('data/modularCapacitiesAndRoutingDualDecomp_lowerBound.txt');
bestLowerBound = load ('data/modularCapacitiesAndRoutingDualDecomp_bestLowerBound.txt');
pie = load ('data/modularCapacitiesAndRoutingDualDecomp_pie.txt');
heuristicLBCost = load ('data/modularCapacitiesAndRoutingDualDecomp_heuristicLBCost.txt');

objFunc
lowerBound

multiplot ('Time', 'Cost' , [objFunc 210*ones(size(objFunc,1),1) lowerBound(:,2)] , [] , [0 250] , '-' , '\pi_e' , pie , [] , [0 2] , '-');

multiplot ('Time', 'Cost' , [objFunc 210*ones(size(objFunc,1),1) bestLowerBound(:,2)] , [] , [0 250] , '-' , '\pi_e' , pie , [] , [0 2] , '-');

bestLowerBound
heuristicLBCost
bestObjFunc(end)

